home
***
CD-ROM
|
disk
|
FTP
|
other
***
search
/
C/C++ Users Group Library 1996 July
/
C-C++ Users Group Library July 1996.iso
/
vol_200
/
245_01
/
lca31.c
< prev
next >
Wrap
Text File
|
1987-10-26
|
16KB
|
576 lines
/* (3,1) Linear Cellular Automaton */
/* Reference: */
/* */
/* Kenneth E. Perry */
/* Abstract Mathematical Art */
/* BYTE */
/* December, 1986 */
/* pages 181-192 */
/* Copyright (C) 1987 */
/* Harold V. McIntosh */
/* Gerardo Cisneros S. */
/* G. Cisneros, 4.3.87 */
/* 10 April 1987 - modified for (4,2) [HVM] */
/* 26 April 1987 - Multiple menus [HVM] */
/* 28 April 1987 - back modified to (4,1) [HVM] */
/* 28 April 1987 - version for XVI Feria de Puebla [HVM] */
/* 14 May 1987 - modified for (3,1) and general rule [HVM] */
# include <bdos.h>
# define COLGRAF 4 /* graph resolution */
# define T80X25 3 /* text resolution */
# define WHCYMAG 1 /* color quad for normal screen */
# define AL 320 /* array length (screen width) */
# define TS 7 /* distinct sums w/totalistic rule */
# define DS 27 /* number of distinct sums */
# define KK 3 /* number of states per cell */
# define NX 21 /* number of sample rules */
char xrule[NX][KK][KK][KK];
char ixrule[NX][DS]=
"001200111011011222201121100", /* interfaces of 2 vel */
"001200211221001010221220120",
"001210120212011110020110100", /* mottled background */
"002110012010121101201012211", /* (0*102211200)* still */
"002110012110021201001012221", /* symmetrized */
"002200011112011011012022200", /* threads on triangles */
"010002101012100210102011211", /* dendrites */
"010020202121101100222101120", /* gldrs, puffr trn */
"010022101200112021010000110", /* threads among trngs */
"010102020102020202020202020",
"010102022102022221022221210", /* open triangles */
"011200122100121212211112212",
"012100200110011022102202022", /* 2 comes and goes */
"012112212012012112012212210", /* 12 irresistible gldr */
"020101011101011112011112120", /* diagonal black gaps */
"020122100021210222210022111", /* black triangles */
"020201010201010102010102020", /* class iv, sparse */
"021111222210112122112122220", /* 102 ia a barrier */
"021211110211110102110102020", /* blue bgrnd cl iv */
"210021220202112222100220222", /* 0122 barrier */
"221212021202011210012120111" /* 2 glider on 1 bkgrnd */
;
char xx[4], rule[DS+1], ascrule[KK][KK][KK];
int binrule[KK][KK][KK], arule[DS], arr1[AL], arr2[AL];
char trule[TS]="0000000";
main()
{
int i, j, i0, i1, i2, jj, n;
int more = 'r';
char a, b, c;
for (i=0; i<NX; i++) { /* copy to 3-index array */
i0=0; i1=0; i2=0;
for (j=0; j<DS; j++) {
xrule[i][i0][i1][i2]=ixrule[i][j];
i2++;
if (i2>2) {i2=0; i1++;};
if (i1>2) {i1=0; i0++;};
if (i0>2) {i2=0; };
};};
videopalette(WHCYMAG); /* white/cyan/magenta */
tuto();
while (!kbdst()) jj=rand(); /* wait for keypress */
c=kbdin(); /* ignore it */
jj=rand()%NX;
for (i0=0; i0<KK; i0++) {
for (i1=0; i1<KK; i1++) {
for (i2=0; i2<KK; i2++) {
ascrule[i0][i1][i2] = xrule[jj][i0][i1][i2]; /* random sample rule */
};};};
rule[DS]=0;
ranlin(); /* random initial array */
videomode(T80X25);
videoscroll(3,0,5,71,0,3); /* menu on blue background */
videoscroll(19,0,24,71,0,3);
xmenu(jj+1); /* show initial rule */
while (more!='n') { /* execute multiple runs */
rmenu();
lmenu();
while (0<1) { /* set up one run */
c=kbdin();
if (c=='g') break; /* go draw graph */
if (c=='q') more='n'; /* quit for good */
if (more=='n') break;
switch (c) {
case '@': /* numbered tot rule */
nutoto(numin(0));
totoasc();
rmenu();
videocursor(0,4,0);
break;
case '$': /* dozen totalistics */
j=numin(0);
for (i=0; i<12; i++) {
nutoto(j+i);
totoasc();
ranlin();
evolve(rule);
};
videomode(T80X25);
rmenu();
lmenu();
break;
case 't': /* totalistic rule */
xblnk();
tmenu();
edtrule();
totoasc();
for (i0=0; i0<3; i0++) {
for (i1=0; i1<3; i1++) {
for (i2=0; i2<3; i2++) {
ascrule[i0][i1][i2]=trule[i0+i1+i2];
};};};
videocursor(0,4,0);
rmenu();
xmenu(totonu(0));
break;
case 'r': /* edit rule */
xblnk();
edrule();
videocursor(0,4,0);
rmenu();
break;
case 'l': /* edit cell string */
xblnk();
edline();
videocursor(0,3,0);
lmenu();
break;
case '#': /* read stored rule */
xmenu(NX);
n=(i=lim(1,numin(0),NX))-1;
xmenu(i);
xtoasc(n);
rmenu();
break;
case 'D': /* run through samples */
for (i=0; i<NX; i++) {
xmenu(i);
xtoasc(i);
ranlin();
evolve(rule);
};
videomode(T80X25);
rmenu();
break;
case 'u': /* sparse init arry */
xblnk();
for (i=0; i<AL; i++) arr1[i]=0;
arr1[AL/4]=1;
arr1[AL/2]=2;
arr1[(3*AL)/4]=1;
arr1[(3*AL)/4+2]=2;
lmenu();
break;
case 'x': /* random rule */
xblnk();
for (i0=0; i0<KK; i0++) {
for (i1=0; i1<KK; i1++) {
for (i2=0; i2<KK; i2++) {
if ((KK*(KK*i0+i1)+i2)%4 == 0) i=rand();
ascrule[i0][i1][i2]='0'+i%3;
i/=3;
};};};
rmenu();
break;
case 'y': /* random line */
xblnk();
ranlin();
lmenu();
break;
case 'Y': /* symmetrize rule */
for (i0=0; i0<KK; i0++) {
for (i1=0; i1<KK; i1++) {
for (i2=0; i2<KK; i2++) {
ascrule[i2][i1][i0]=ascrule[i0][i1][i2];
};};};
break;
case 'B': /* begin barrier */
a=kbdin();
b=kbdin();
ascrule[0][a-'0'][b-'0']=a;
ascrule[1][a-'0'][b-'0']=a;
ascrule[2][a-'0'][b-'0']=a;
rmenu();
break;
case 'E': /* end barrier */
a=kbdin();
b=kbdin();
ascrule[a-'0'][b-'0'][0]=b;
ascrule[a-'0'][b-'0'][1]=b;
ascrule[a-'0'][b-'0'][2]=b;
rmenu();
break;
case 'L': /* left glider link */
a=kbdin();
b=kbdin();
c=kbdin();
ascrule[a-'0'][b-'0'][c-'0']=c;
rmenu();
break;
case 'R': /* left glider link */
a=kbdin();
b=kbdin();
c=kbdin();
ascrule[a-'0'][b-'0'][c-'0']=a;
rmenu();
break;
case 'S': /* still life link */
a=kbdin();
b=kbdin();
c=kbdin();
ascrule[a-'0'][b-'0'][c-'0']=b;
rmenu();
break;
case '=':
for (i=1; i<8; i++) {
for (j=0; j<40; j++) arr1[40*i+j]=arr1[j];};
lmenu();
break;
case '~':
for (i=1; i<16; i++) {
for (j=0; j<20; j++) arr1[20*i+j]=arr1[j];};
lmenu();
break;
default: break;
};
};
if (more=='n') break;
do {
evolve(rule);
videocursor(0,0,0);
scrstr("More?");
videocursor(0,0,34);
scrstr("y/n/cr");
more=kbdin();
} while (more=='\015');
videomode(T80X25); /* reset the screen */
if (more=='n') break;
};
videomode(T80X25);}
/* edit the rule */
edrule() {
char c;
int i, i0, i1, i2;
i=6; i0=0; i1=0, i2=0;
while (0<1) {
videocursor(0,3,i);
c = kbdin();
if (c == '\015') break; /* carriage return exits */
switch (c) {
case '0': case '1': case '2': /* state */
ascrule[i0][i1][i2] = c;
i2++;
if (i2>2) {i2=0; i1++;};
if (i1>2) {i1=0; i0++;};
if (i0>2) {i2=0; };
videocattr(0,c,3,1);
if (i<6+DS) i++;
break;
case ' ': case '\315': /* space = advance */
i2++;
if (i2>2) {i2=0; i1++;};
if (i1>2) {i1=0; i0++;};
if (i0>2) {i2=0; };
if (i<6+DS) i++;
break;
case '\010': case '\313': /* backspace */
if (i2!=0) i2--; else {i2=2;
if (i1!=0) i1--; else {i1=2;
if (i0!=0) i0--; else {i0=2;
};};};
if (i>6) i--;